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ABSTRACT 

We present results from our survey for planetary transits in the field of the intermediate 
age (~2.5 Gyr), metal-rich ([Fe/H]~+0.07) open cluster NGC 6819. We have obtained high- 
precision time-series photometry for over 38,000 stars in this field and have developed an 
effective matched-filter algorithm to search for photometric transits. This algorithm identified 
8 candidate stars showing multiple transit-like events, plus 3 stars with single eclipses. On 
closer inspection, while most are shown to be low mass stellar binaries, some of these events 
could be due to brown dwarf companions. The data for one of the single-transit candidates 
indicates a minimum radius for the companion similar to that of HD 209458b. 

Key words: methods: data analysis - open clusters and associations: general - open clusters 
and associations: individual: NGC 6819 - binaries: eclipsing - stars: low mass, brown dwarfs 
- planetary systems. 



1 INTRODUCTION 

Following the discoveries of Hot Jupiter-type extra-solar planets, 
it was apparent that a photometric survey to detect these plan- 
ets in transit was feasible from ground-based telescopes. This 
was confirmed by the de tecti on of transits by HD 209458b by 
ICharbonneau et all fcOOOl) and lHenrv et alj feOOCt) . Such a survey 
has the potential to broaden the search for planets over a greater 
volume of space by monitoring stars at fainter magnitudes than are 
accessible to the spectroscopic radial velocity technique. The tran- 
sit method enables us to probe far greater numbers of stars simulta- 
neously and so obtain a statistically significant sample of planets in 
a comparatively short time. Discoveries can be used to determine 
the abundance of planets in a range of stellar environments, allow- 
ing us to investigate the relationship between planet formation and 
key properties; for example metallicity, age, stellar and radiation 
density. Follow-up observations of transiting planets can provide 



* Based on observations made with the Isaac Newton Telescope operated 
on the island of La Palma by the Isaac Newton Group in the Spanish Ob- 
servatorio del Roque de los Muchachos of the Instituto de Astrofisica de 
Canarias. 



vital information on individual planets - their true masses, radii 
and orbital inclinations - crucial for testing theories of planetary 
structure. Once transit candidates are obtained from photometric 
surveys, these data combined with spectrographic observations will 
not only provide this information but will also distinguish planetary 
transits from grazing-incidence eclipses by stellar and brown dwarf 
companions. Follow-up observations are discussed in greater detail 
in Section lTTI 



From the radial velocity survey results, we know that approx- 
imately 1 percent of F-K-type solar neighbourhood stars harbour 
hot Jupiters (giant planets with periods of ~ 2-6 days). The typi- 
cal orbital separation (0.05 AU) for these planets implies that about 
10% of them should exhibit transits if their orbital planes are ran- 
domly oriented with respect to the line of sight. Assuming that 
planetary orbits are distributed isotropically, we can expect there to 
be one transit for every 1000 stars. The detection probability is af- 
fected by factors such as the sampling rate of the observations, the 
efficiency of detection from the data and the fraction of late-type 
dwarfs in the sample. Correspondingly, transit surveys that moni- 
tor tens of thousands of stars simultaneously may expect to detect 
tens of planets. This requirement must then be balanced against the 
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Table 1. Basic data on NGC 6819. Data taken from Kalirai et al. (2001) and 
the SIMBAD database. 



RA (J2000.0) 


19 h 41 m 


Dec (J2000.0) 


+40° 11' 


I 


73°. 97 


b 


+8°. 48 


Distance (pc) 


2754 ± 305 


Radius 


~9.5' 


Age (Gyr) 


2.5 


[Fe/H] 


+0.07 


E(B-V) 


0.10 



need to avoid fields so crowded that blending makes stars difficult 
to m easure precis ely. 

Janes ( 1996) suggested that open clusters would make a good 
compromise for ground-based surveys. They provide large numbers 
of stars within a relatively small field with minimal blending. They 
also provide a distinct population of stars of known age and metal- 
licity and allow us to probe stars in the environment where they 
formed. Additionally, these fields provide a separate population of 
background stars for comparison. 

A number of groups are pursuing transit surveys using large 
(2^1m) telescopes with wide field, mosaic CCD cameras. Most no- 
tably, the EXPLORE project has used the CTIO-4m and the 3.6m 
CFHT to observe two Galactic plane fields in 2001 lYee et alJ 
2002). These data have revealed 3 possible planetary transit can- 
didates, and the t eam are in the proces s of obtaining radial velocity 
follow-up ( Malle ~Ornelas et alj2002l) . The OGLE group have used 
their microlensing observations of Galactic disc stars to search for 
transits in the lightcurves o f ~52,000 star s , yield ing 59 candidates 
so far lUdalski et alj2002h . lDreizler et alJ 12002) obtained classifi- 
cation spectra for 16 of these stars, allowing them to estimate the 
radius of the primary and infer the radius of the companion. This 
analysis ruled out 14 candidates as having stellar-mass compan- 
ions, while the companions of two objects were found to have radi i 
similar to that of HD 209458b. Recently. iMocheiska et ail 1 2002) 
have undertaken a survey of open clusters using the F. L. Whipple 
Observatory's 1.2m telescope, discovering 47 new, low amplitude 
variables. 

In 1999 we began a survey of three open clusters for planetary 
transits: NGC 6819, 7789 and 6940. The wide field of view of the 
Isaac Newton Telescope's (INT) Wide Field Camera (WFC) is ide- 
ally suited to this task, and we were awarded a total of 3 bright runs 
of 10 nights each in June/July 1999 and September 2000 for these 
observations. 

Open cluster NGC 6819 was observed during the first 19 of 
these nights and the results presented here stem from our analysis of 
these data. The basic parameters of this cluster are given in Tabled 
Previo us relevant work on this cluster was dis cussed inlStreet et alJ 
(2002); most important is a recent study bv lKaliraTeUriTT2fj(j1I) 
which provided B and V magnitudes for large numbers of stars in 
this field. 

In the rest of this paper, Section|2|details these observations, 
and Sections [5] and [4] present the data reduction procedure. Sec- 
tion|5|discusses our transit-detection algorithm while our results are 
presented in Section|6| Finally we draw conclusions in Section[8| 



Table 2. Gain and readout noise values for the four WFC CCDs, taken from 
the Cambridge Astronomical Survey Unit webpage. 



CCD No. 


Gain(e" ADLT 1 ) 


Readout noise (ADU) 


1 


3.12 


7.9 


2 


3.19 


6.4 


3 


2.96 


8.3 


4 


2.22 


8.3 



2 OBSERVATIONS 

We observed NGC 6819 on 19 nights during 1999 June 22-30 
and 1999 July 22-31, using the 2.5m Isaac Newton Telescope, La 
Palma. The Wide Field Camera (WFC) employs four 2048x4096 
pixel EEV CCDs to image a ~0.5° x0.5° field of view with a pixel 
scale of 0.33" pix" 1 . The gain and readout noise values for each 
CCD were taken from the Cambridge Astronomical Survey Unit 
webpage 1 and are listed in Table|2| 

Three open clusters, NGC 6819, 6940 and 7789, were ob- 
served in rotation during these two runs, taking pairs of 300s ex- 
posures through the Sloan r' filter on each visit. The readout time 
of the WFC at the time required 160 sec of deadtime between ex- 
posures. No dithering was applied between exposures as we aimed 
to place each star as close to the same pixel each time as possi- 
ble. In practise, the x, y shifts between images were up to a few 
pixels. NGC 6819 was observed for ~ 7 — 8 hours each night, typ- 
ically resulting in 16 - 25 frames per night or about 2 frames/hour. 
In total, 361, 384, 364 and 325 frames were obtained of this clus- 
ter for CCDs 1-4 respectively. The number of available frames 
varied due to unpredictable readout failures which could affect in- 
dividual CCDs. The average gap between pairs of exposures was, at 
most, roughly an hour and we had good observing conditions on all 
nights. The field of view covered by the WFC is shown in FigureQ 



3 DATA REDUCTION 

We developed a semi-automated data reduction pipeline in order 
to proces s this large dataset. This pipeline h as been previo usly de- 
scribed in Street et al. 120021) and in detail in IS treed 120021) . 

The debiasing, flatfielding and trimming of the fram es were 
carried out using the Staiiink package FIGARO Shortrid ge et all 
1 1998). A correction for the nonlinear response of each chip was 
also applied. Point-spread function (PSF) p hotometry was per- 
formed using IRAF's DAOPHOT task lStetson ll987l) . It was found 
that stars subtracted from an image using a fixed PSF showed resid- 
uals that varied with position. These residuals were best reduced by 
employing DAOPHOT's "penny 2" function and allowing it to vary 
quadratically with position. This is a two-component model, con- 
sisting of an elliptical Gaussian core and Lorentzian wings. Both 
parts of the model are aligned along separate and arbitrary posi- 
tion angles. DAOPHOT handles pixel defects, cosmic rays etc. by 
employing a formula which reduces the weights of pixels that do 
not converge towards the model as the fit is calculated. This is dis- 
cussed in more detail in IPavisI ll994l) . The post-processing (dis- 
cussed below) is also able to detect and remove strongly outly- 
ing points. Where a star of interest is found to lie close to dead 

1 http://www.ast.cam.ac.uk/~wfcsur/ccd.html 
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Figure 1. The WFC ~0.5 x0.5 field of view covering open cluster NGC 6819; North is to the top of the image while East is to the left. The CCD numbers 
are given in the corners. 



columns/pixels conclusions have been drawn with caution. We 
chose to have the star positions re-fitted independently in each 
frame, having found that the cross-correlation technique aligned 
the star centroids to around ~ 1 pixel accuracy. 

For a significant number of images, we found that a position- 
dependent element still remained in the magnitude residuals, partic- 
ularly dominant along the long (y) axis of the CCD. To counteract 
this problem, we have employed our own post-processing software, 
described in lStreet et alj 120021) . This included a procedure which 
cross-correlates all lightcurves in the sample in order to identify 
and remove remaining systematic trends. 

Following this processing, the precision achieved is illustrated 
by plotting the RMS scatter in each star's lightcurve against its 
weighted mean magnitude over the whole dataset. Figure|2|shows 
these plots for each CCD, and for reference shows the effects of 
the main expected sources of noise. While some systematic effects 
remain in the data, our software improves the precision particu- 
larly at brighter magnitudes, where we can achieve the ~0.004 mag 
precision required to detect planetary transits. We notice that the 
residual systematic variations are reduced to a level of ~0.0025™ 
0.0035 mags. After post-processing, these residuals do not appear 
to show a positional distribution. We also notice that the "back- 
bone" of points falls slightly below the theoretical noise prediction 
at the faint end. This seems to be due to DAOPHOT underestimat- 
ing the magnitude errors used to weight the calculations at fainter 
magnitudes. We note the presence of "clumps" of stars with high 
RMS in CCD3. Investigation of these points reveals that they are 
hot pixels located in the vignetted areas and dead columns which 
this CCD suffers from. 

Astrometric positions for all the stars i n our samp le were ob- 
tained using the method described in IStreet et aljfe002l) . The aver- 
age RMS error in the resulting RA and Dec are presented in Table|5| 
and correspond to an RMS scatter of less than 1 pixel on the CCD. 



Table 3. Average RMS error on RA and Dec positions in arcsec for each 
CCD. 



CCD No. 


<5RA 


<5Dec 




(arcsec) 


(arcsec) 


1 


0.118 


0.110 


2 


0.232 


0.363 


3 


0.289 


0.294 


4 


0.172 


0.180 



4 STELLAR RADII 

4.1 Colour Index Calibration 

iKalirail 1200 ll) kindly provided us with B — V colours indices 
for many of the stars in the NGC 6819 field. These data and the 
proced ure used to cross-identify stars are described in lStreet et alJ 
(2002). Although most of our data were taken in the Sloan r band, 
we obtained enough Sloan i data to calibrate approximate V — R 
colours for all the stars in our sample in the following way. We cal- 
culated the mean broadband flux using known passband functions 
and the Bruzual, Persson, Gunn and Stryker atlas of stellar spec- 
tra. The mean broadband flux was then used to calculate theoretical 
magnitudes and colours for a range of stellar spectral types relative 
to the flux from Vega in that passband. The INT instrumental Sloan 
r — i colours were calibrated by superimposing the XCAL plot of 
Sloan r — i against B — V over that of the INT stars, and applying 
vertical and horizontal offsets. The horizontal offset provided the 
calibration factor for the INT data, in the sense that the true Sloan 
r — i colour of a star (r — i) is found from the instrumental one 
(r — i)inst by adding the offset, A, — 

r - i = (r - i)mst + A r -i. (1) 
A r _ t was found to be (T50, (T50, (T30 and (T41 for the 
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Figure 2. The variation of RMS scatter in star lightcurves with magnitude for each of the four CCDs. Superimposed curves show the effects of various sources 
of noise. 



data from CCDs 1, 2, 3 and 4 respectively. In all cases, the vertical 
offset (the difference between the theoretical and measured B — V) 
was found to be ffi"2. This we attribute to extinction in the direction 
of the cluster and is not very different from the value of E{B — 
V) = 0.1 measured bv lKalirai et aflbOOll) . 

We used a similar method to convert the now-calibrated Sloan 
r — i colours into Johnson V — R. XCAL was used to produce a 
dataset of V — R and corresponding Sloan r — i values. To derive 
a formula to convert Sloan r — i colours into Johnson V — R, a 
function was fitted to these data using the method of least squares. 
As the shape of the curve changes at r — i ~ 0.4, two functions 
were fitted; a straight line for r — i values ~ —0.2 ~ 0.4 and an 
exponential function for the remaining curve: 



4.2 Colour-Magnitude Diagram 



V-R = 



0.009 + 0.93(r -i) 
1.30 - 1.72e- 1 - 6(r - <) 



-0.2 < r - 
0.356 < r 



i < 0.356 
-i < 2.5 



(2) 



These relations were then used to calculate V 
dices for all INT stars with Sloan r — i colours. 



R colour in- 



Figure|5|shows the (V, V — R) colour-magnitude diagrams (CMD) 
for the four CCDs. The cluster main sequence is clearly visible in 
the data from CCD4 and faintly in the other three plots. This is ex- 
pected from the radius of the cluster Kalirai et alJ 1200 ll) which fits 
within the field of view of one of the WFC CCDs. Field stars greatly 
outnumber cluster members in these data, and are located above 
and below the cluster main sequence in the CMD. The field stars 
above and below the cluster main sequence are predominantly main 
sequence stars at closer and more remote distances than the cluster, 
respectively. In our following analysis, we assume that all stars in 
our sample are main sequence. This is reasonable since class IV 
subgiants are sufficiently rare that their frequency in the sample is 
negligible. Giant stars are so bright (-0^4 < M v < 1™2) that 
for one to be measurable (unsaturated) in our data it would have to 
be at a distance of >9Kpc. We therefore adopt main sequence re- 
lationships and use the likely distance and spectral type and hence 
radius of the stars in our sample. 
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Single transit candidates 
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(a) CCD 1 
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(b) CCD 2 



CCD3 data 




0.2 0.4 0.6 0.8 1 1.2 1.4 0.2 0.4 0.6 0.8 1 1.2 1.4 

V-R V-R 

(c) CCD 3 (d) CCD 4 

Figure 3. V.vs.V — R colour- magnitude diagrams for the four WFC CCDs. The transit candidate stars are highlighted. 
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Figure 4. RMS scatter against V — R colour for all ^38,000 stars in our sample. The overlaid curves give the predicted transit depth for cluster member stars 
transited by planets of various radii. 
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4.3 Stellar Radii 

We use the V — R index to estimate the radius of each star in this 
survey. This was done by interpolating between measured values 
of V — R a nd R t for the range of main sequence star types from 
lGravHl992t) . These values were s upplemented at the low-mass end 
by data from lReid & Gizisl < 19971) . who list V — R index, absolute 
magnitude and spectral type for 106 low-mass systems. The abso- 
lute magnitudes were then used to calculate the radii of these stars. 
An exponential curve was least-squares fitted to the Reid and Gizis 
data and used to calculate values of stellar radii at fixed intervals of 
V — R in order to provide one smooth, continuous dataset. Inter- 
polation over this dataset was then used to compute main sequence 
star radii from the V — R colour index of each star. Figure[4]plots 
RMS scatter versus V — R colour. This is overlaid with curves il- 
lustrating the predicted transit depth of planets with radii of 0.5, 1.0 
and 2.0 Rj up orbiting cluster member stars of various masses. If the 
RMS of a given star falls below one of these curves, then we would 
expect to be able to detect a transit of a planet that size around that 
star. 

We find that ~30% of the stars in our sample fall below the 
1.0 Kj U p -transit curve, while ~79% fall below the 2.0 Rj up line. 
Thus around ~ 11,500 and ~30,000 stars, respectively, are mea- 
sured to sufficient precision to allow the detection of transits. As- 
suming ~1% of all main sequence stars have hot Jupiter compan- 
ions and ~10% of those transit, then we can roughly expect to de- 
tect ~1 1 Jupiter-radius objects in our data. 



5 TRANSIT DETECTION ALGORITHM 

Having obtained the required high-precision photometry, there are 
a number of different approaches to the problem of detecting transit 
events. We have developed our own transit-finding software, using 
the method of matched-filter analysis. 

After identifying and removing known large-amplitude vari- 
ables from the data, the software works in two stages. The first 
stage, or "standard search", generates a series of model lightcurves 
with a single transit. These models are generated for a range of 
transit durations 0.5 < St < 5.0 hours, in intervals of 0.25 hours 
and with the time of mid-transit ranging from the start to the end 
of the observing campaign in steps of St/ 4. A constant magnitude 
is least-squares-fitted to each lightcurve and the corresponding Xc 
is calculated. Each model is then \ 2 -fitted to each lightcurve, the 
transit depth and out-of-transit magnitude being optimised by min- 
imising the quality-of-fit statistic Xm- The details of the best-fitting 
model with the lowest value of Xm are stored. 

A transit-finder index, Ax?/, is then calculated as: 

Ax?/=Xc-Xm- (3) 

Plotting this index against Xm allows us to separate transit 
events from constant stars and other types of variables. A genuine 
transit event would be expected to show a significant improvement 
in x 2 when comparing the fit of a constant line and a suitable transit 
model; hence it would have a relatively low value of Xm and a 
high value of Ax?/- To isolate these candidates, a straight line is 
least-squares fitted to the "backbone" of points. This fit is iterated, 
rejecting all points ±3u above the line, until the parameters change 
by less than 0.0001. A cut-offline is established by raising this line 
by +Na where TV is set by the user; all stars falling to the top- 
left of this cut-off are regarded as candidates. To illustrate how this 



isolates transit candidates, we tested the algorithm by injecting fake 
transits with known parameters into the data stream. Transits were 
added to 1% of stars in the CCD1 data, with a period of 3.4days, 
a duration of 2.5 hours, and an amplitude of I "02. The first transit 
occurred at HJD 2451355.5 and at every multiple of the period. The 
transit-finder algorithm was then applied to this modified dataset 
and the resulting plot of Ax?/ against Xm is shown in Figure [5] 
The transits of an HD 20945 8b-like planet are clearly separated 
from the rest of the data. This figure was used to set the detection 
threshold. It was found that a +4<j threshold retains all but 5 of the 
90 fake transits while excluding ~99.2% of the constant stars. 

The second stage of the search is a "period search" applied to 
candidates highlighted by the standard search. For each candidate, 
multiple-transit models are generated across a range of periods (2- 
5 days) and fitted to the lightcurve as described above. Once again, 
the minimum value of Xm sets the best-fitting model, and candi- 
dates are selected by the method described above. This two stage 
approach ensures that all the relevant transit-parameters are deter- 
mined for all candidates, but restricts the number of least-squares 
fits required by applying the period search to transit-candidates 
only. This allows a statistically-optimal matched-filter technique to 
be applied without prohibitive computational time requirements. 

We note that an element of human judgement enters into our 
transit detection procedure. The candidates presented by the algo- 
rithm are sorted by manual examination. In the process we have 
rejected a large number of "possible" transits; lightcurves which 
show dips at the beginning or end of a night which do not repeat, 
for example, or those that show dips sampled with very few data- 
points. Of course this means we could potentially miss transit in- 
gresses/egresses, but a real candidate must show at least two well- 
sampled transit events. 



6 RESULTS 

In total, over 38,000 star lightcurves have been analysed in this way. 
The transit search algorithm highlighted 276 stars worthy of fur- 
ther investigation, and these were examined manually. The major- 
ity (~51.9%) were found to show only a few, fainter-than-average, 
scattered points. The cause of the scattering was found to be one 
of three situations: (a) the presence of nearby or blended com- 
panion star(s), (b) the star is bright and saturated in a significant 
number of images or (c) the star falls close to a dead column or 
vignetted area on the CCD. The algorithm also detected what we 
judge to be stellar variability in ~20.3% of cases; most of these 
stars showed eclipses due to stellar companions while some dis- 
played low-amplitude "dips" in brightness due to stellar activity, 
although on longer timescales than transits. No obvious explana- 
tion for spurious detection could be found for 67 of these stars; 
in these cases examination of the lightcurve revealed unconvinc- 
ing "transits" consisting of well-scattered points, often on nights of 
poor conditions. 

Of the remaining stars (~2.9%), 8 appear to show short- 
duration transit-like eclipses. The sample includes a number of ac- 
tive stars which show brief eclipses. This is not unexpected since 
transit amplitude scales inversely with star radius squared, while 
stellar activity is more common among small, young stars. All 
short-duration eclipses were considered regardless of amplitude, 
since hot Jupiter transits could reach depths of up to several tenths 
of a magnitude, given a late-M type primary. 

A full lightcurve solution to a fitted model is not possible be- 
cause of the sparseness of the data. However, the V — R colours 
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Figure 5. Plot of Ax^ against \m f° r me CCD1 dataset with fake HD 209458-like transits added; these points are highlighted. Transits fall to the upper left 
of this diagram, above the main backbone of points - the threshold is set at +4cr. 
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Figure 6. Plots of t against \m f° r a " f° ur CCDs. All lightcurves above the superimposed cut-off thresholds are visually examined for transits. 
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for these stars can be used to estimate the radius of the primary star 
(R,), assuming for the moment that the star is main sequence and 
that negligible light is contributed from the companion body. The 
amplitude, 8m, of the transit is proportional to the ratio of the star's 
radius to that of the companion (R c )'- 

We estimate R c from Equation|4| noting that this gives a lower 
limit because a larger companion can cover the same fraction of pri- 
mary star if the eclipse is partial rather than total. The radius of the 
companion gives a general indication of its nature. However, while 
the radius of a main sequence M star can be ~0. 1-0.5 Rq , the radii 
of gas giant planets (R J tlp ~0. 1 Rq ) are thought to be similar to 
those of brown dwarfs (<0.2 R0 ), owing to their degenerate na- 
ture. The transits of HD 209458b have also shown that hot Jupiter 
radii can be larger than expected due to their early proximity to their 
primary star slowing the rate of contraction iBurrows et al] §000). 
For this reason, follow-up radial velocity measurements yielding 
the minimum mass will be required to distinguish between plane- 
tary and stellar companions. 

Table [4] presents the details of the 8 candidates. Note that 
one candidate is marked as being blended with nearby stars. The 
conclusions drawn from these stars come with the caveat that 
the results need to be confirmed. All the candidates are found to 
have a minimum companion radius below 0.5 Rq , while 3 have 
R c < 0.25 Rq . The phased lightcurves of the stars with M-dwarf 
or smaller companions are presented in Figure The candidates 
are discussed individually below. Some of the companion objects 
could be brown dwarfs although most are found to be low mass 
stars. In these cases, we note that stellar companions will con- 
tribute by reddening the measured colour - this would mean that 
the primary and companion stars are of larger radii than calculated 
here, as would inclinations of less than 90°. Stellar binaries would 
also exhibit secondary eclipses not seen in transit lightcurves. We 
have not included objects which clearly show eclipses of differ- 
ent depths, as these will be stellar binaries. Where objects showing 
similar eclipse-depths turn out to be stars, then the period will be 
twice that given in Table|4| 

6.1 Star 249 - P=2.233 days, 5m=(T19 

With all binary objects with periods as short as these it is possible 
for the rotation of the primary to have become synchronised with 
the orbital period of the companion. If this is the case, then stellar 
activity on the primary is expected, resulting in the variable out-of- 
eclipse lightcurves. This seems to be the case for star 249: the im- 
plied companion radius (0.37 R0 ) is that of an M-dwarf or larger. 
While the eclipses are not well sampled there is some suggestion 
that they may be rounded-bottomed. 

6.2 Star 4619 -P=3.682 days, <5m=(T03 

This candidate shows the classic transit lightcurve: sharp 
ingress/egress to low amplitude eclipse with no out of transit vari- 
ations. The eclipse profiles are not well sampled but could be 
rounded-bottomed, and the period is typical of the known hot 
Jupiters. The V — R colour (0^283) indicates a primary radius of 
~1.5 R© which together with the low amplitude implies a compan- 
ion radius of ~0.26 R0 . The companion could be a brown dwarf. 



6.3 Star 6690 - P=1.682 days, <5m=<T09 

This candidate also shows the expect transit lightcurve except that 
the eclipses appear to have a sharp, pointed profile, suggesting that 
these are grazing incidence eclipses. This would mean that the com- 
panion radius is larger than 0.23 Rq . However, the colour (0^719) 
and amplitude (0^09) imply a relatively small primary and sec- 
ondary radii. In this case the companion is likely to be a low mass 
star. 



6.4 Star 10400 - P=1.46 days, 6m=(f? 14 

Measurements of this star are complicated by the presence of close, 
blended companions. The colour (0^343) and amplitude (0YM4) 
indicate that both primary and secondary radii are stellar. 



6.5 Star 1 1644 - P=2.302 days, 5m=tt? 04 

The lightcurve of this star also shows some modulation between 
the eclipses, which seem to be flat-bottomed; further photometric 
data are needed to confirm this. The colour (0^397) and amplitude 
(0^04) suggest that the primary has a radius of ~1.32R0 while the 
minimum secondary radius is found to be ~0.264 Rq . The com- 
panion object could therefore be a brown dwarf. 

6.6 Star 16155 - P=3.486 days, < 5m=0 m 07 

This lightcurve is similar to that of star 6690, with pointed eclipse 
profiles. However, in this case even the companion's minimum ra- 
dius (0.35 Rq ) implies a small star and a grazing incidence eclipse 
suggests a larger companion. 



6.7 Star 20910 - P=1.3112 days, 5m=(Tl 

This lightcurve displays eclipses apparently rounded and weak ro- 
tational modulation of the out-of-eclipse lightcurve. The 1.3 day 
period and early-K spectral type suggest that magnetic starspot ac- 
tivity driven by the tidally-synchronised rotation of the primary 
is responsible for the modulation. The lack of a clear secondary 
eclipse suggests a very low effective temperature for the compan- 
ion, which is probably a late-M dwarf or a brown dwarf. 

6.8 Star 22790 -P=3.621 days, <5m=0 m 25 

The eclipse profile again suggests a grazing incidence orbit while 
the minimum radius (0.37 Rq ) implies that the companion is a low 
mass star. 



6.9 Single-transit Candidates 

Normally we would require at least two transits in a lightcurve in 
order to consider a star as a candidate, but the range of possible 
orbital periods means that in ~50% of cases, a transiting planet 
will only show a single transit in 20 nights of observations. The 
procedure outlined above also identified 3 lightcurves which appear 
to show single transit-like eclipses. The calculated minimum radii 
of these companions are all less than 0.3 Rq . Table [5] gives the 
details of these stars, while Figure [8] displays the full lightcurves 
next to lightcurves of the "transits". 
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Table 4. The parameters of the stars which show transit-like eclipses. 



Star 


V 


V -R 


8m 


St 


R* 


Rc 


Period 


Epoch 


N t r 


RA 


Dec 




(mag) 


(mag) 


(mag) 


(hours) 


(R©) 


(R©) 


(days) 


(HJD-2400000) 




(J2000.0) 


(J2000.0) 


249 


18.879 


0.623 


0.19 


2.4 


0.84 


0.37 


2.233(31) 


51352.535(1) 


2 


19 42 15.05 


+40 04 42.1 


4619 


16.603 


0.283 


0.03 


4.8 


1.51 


0.26 


3.682(1) 


51387.609(4) 


2 


19 41 21.32 


+40 02 14.3 


6690 


18.667 


0.719 


0.09 


3.1 


0.77 


0.23 


1.682(1) 


51356.578(7) 


3 


19 40 56.71 


+40 05 05.0 


10400 a 


18.906 


0.343 


0.14 


4.3 


1.16 


0.43 


1.46(6) 


51357.506(1) 


5 


19 40 05.30 


+40 14 17.7 


11644 


19.130 


0.397 


0.04 


3.6 


1.32 


0.26 


2.302(2) 


51382.534(4) 


3.5 


19 40 13.93 


+40 11 21.9 


16155 


18.018 


0.408 


0.07 


2.6 


1.32 


0.35 


3.486(5) 


51359.526(2) 


3 


19 40 12.44 


+40 00 45.1 


20910 


18.464 


0.798 


0.10 


1.9 


0.73 


0.23 


1.3112(6) 


51356.4832(9) 


5 


19 41 57.10 


+40 18 25.3 


22790 


20.075 


0.813 


0.25 


4.6 


0.73 


0.37 


3.621(2) 


51383.490(15) 


2 


19 41 33.86 


+40 26 35.0 


a Blended 



Table 5. The parameters of the stars which show single transit-like eclipses. 



Star 


V 


V -R 


8m 


St 




Rc 


Epoch 


RA 


Dec 




(mag) 


(mag) 


(mag) 


(hours) 


(R©) 


(R© ) 


(HJD-2400000) 


(J2000.0) 


(J2000.0) 


829 


20.178 


0.500 


0.04 


2.4 


1.08 


0.22 


51385.561(9) 


19 42 07.06 


+39 59 38.9 


8153 a 


21.720 


1.120 


0.21 


4.8 


0.63 


0.27 


51390.522(9) 


19 40 37.98 


+40 01 04.8 


93296 


16.743 


0.534 


0.03 


2.4 


1.00 


0.17 


51389.656(3) 


19 40 21.71 


+40 04 10.0 



a Blended, b Near saturation. 



In all three lightcurves the transit-like events could be flat- 
bottomed, although better sampled photometry is required to de- 
termine this conclusively and to confirm the events. The calculated 
minimum radii suggest brown dwarf companions; the radius of star 
9329 could even be planetary. We strongly urge follow-up of these 
candidates. 



7 FUTURE OBSERVATIONS 

7.1 Follow-up of Transit Candidates 

Since hot Jupiter p lanets can have simil ar radii to brown dwarfs 
and even small stars lBurrows et aflfeOOd) , it is necessary to obtain 
radial velocity measurements in order to determine the minimum 
mass of the companion. Together with high precision, continuously 
sampled lightcurves, the true companion mass can then be derived. 
While this survey has not produced any clear planetary candidates, 
several of our companion objects may be brown dwarfs, and follow- 
up study would be valuable to confirm or deny this. Although the 
candidates from this survey are much fainter than those covered by 
the radial velocity planet hunting surveys, spectroscopic follow-up 
of low-radius companions is possible and can provide very useful 
information. 

Firstly, low resolution spectra of each would provide a more 
secure spectral classification (and radius) of the primary than our 
present estimates based on broad-band colour indices. For the 
faintest candidates (V > 187*5), this will be the only spectroscopy 
follow-up possible. For most of our candidates however, it is possi- 
ble to obtain radial velocity measurements using 8-10m class tele- 
scopes (see for example, lYee et alJ i2002l) '). While they will not be 
precise enough to measure a planetary mass, they will place useful 
limits on the mass, confirming or ruling out stellar companions. 

Continuously sampled photometric follow-up is highly desir- 
able, in two colours if possible. In our original observation strategy 



we decided to cycle around three separate cluster fields in order to 
cover as many stars as possible. In retrospect, we find that continu- 
ous observations of a single field is preferable, in order to get clear, 
well defined eclipses. Highly sampled photometry would clearly 
reveal the eclipse morphology, distinguish total eclipses from graz- 
ing incidence events and allow detailed models to be fitted. For the 
fainter candidates with no radial velocity observations, this will be 
crucial in determining the nature of the system. Photometry would 
also improve the ephemeris, allowing us to time radial velocity ob- 
servations better. The INT/WFC or similarly equiped 2m class tele- 
scope could be used for this purpose. 

7.2 Transit Search Strategy 

The results reported here were obtained from our survey's first 
observing season, and several improvements to our strategy have 
now been adopted as a result. Firstly, well sampled lightcurves are 
crucial. Originally, we tried to include as many stars as possible 
by covering 3 clusters in rotation, resulting in a sampling rate of 
about 2 observations/star/hour. A planetary transit, of typical dura- 
tion ~2.5 hours would be represented by perhaps 4—6 datapoints. 
We have shown that our algorithm can detect such transits (see Fig- 
ure^. However, in practise a greater signal-to-noise is very desir- 
able. It also helps in distinguishing stellar eclipses from transits, 
in determining the properties of the system, and not least in cal- 
culating an accurate ephemeris for follow-up. We have now begun 
continuously sampled observations. 

From the point of view of detecting transits, blended stars in 
crowded fields represent the most significant problem. The addi- 
tional scattering caused in the lightcurve can resemble a transit suf- 
ficiently well to distract the algorithm, and even visual inspection. 
Better sampled data will help to alleviate this, and we are investigat- 
ing image-subtraction techniques which should deal with blending 
more effectively (see for example, iMocheiska et alJ J2002I) ). The 
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Figure 8. The lightcurves of stars showing single transit-like eclipse events. In the left-hand column, arrows mark the location of the suspected transit. 



other major source of false detections is stellar activity and eclips- 
ing binaries, which obviously the algorithm is very good at find- 
ing. We are currently investigating improvements which will reject 
these stars automatically. 

Another issue raised by this work was cluster membership. 



Firstly, it is difficult to know whether any given star is a member 
or not. Although this can be decided by astrometry, few clusters 
have been studied in this way, and usually not to faint enough mag- 
nitudes. The best photometric solution is to obtain good quality 
colour data for colour-magnitude and if possible colour-colour dia- 
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grams. These together with separation from cluster centre measure- 
ments can be used to assign membership probabilities. Secondly, a 
transit survey needs to cover large numbers of stars in its chosen 
population. This survey found that only ~6% of the stars measured 
were cluster members from their colours, amounting to just over 
2113 stars out of 38,118. This total could be improved slightly by 
selecting larger radius clusters which would better cover the field 
of view and reducing the number of unmeasured stars due to blend- 
ing/crowding. A combination of short/long exposures would cover 
stars over a larger range of magnitudes, although some caution is 
required not simply to increase the number of unsuitable early-type 
stars in the sample. Ultimately, however, open clusters only contain 
a few thousand stars at most. This highlights the need to survey a 
number of clusters, in order to observe enough cluster stars with 
similar ages and metallicities to be able to make definitive state- 
ments about the planetary population. The ultimate aim is then to 
extend the survey to include a significant number of clusters cover- 
ing ranges of age and metallicity, which will reveal the dependence 
of planetary formation and evolution on these parameters. As each 
cluster requires around 20 nights of observing time on a 2^-m class 
telescope, these aims are best achieved by collaborative efforts be- 
tween survey teams in order to obtain sufficient telescope time or 
else a large dedicated telescope. 



8 CONCLUSIONS 

We have obtained high-precision photometry on over 38,000 stars 
in the field of open cluster NGC 6819. We have developed an algo- 
rithm which can effectively identify transit-like events in sparsely- 
sampled data. This has produced 8 candidates showing multi- 
ple transit-like events plus a further 3 candidates showing single 
eclipses. Closer analyses of these lightcurves indicates some of 
these candidates could be brown dwarfs, while one has a minimum 
radius similar to that of HD 209458b. Follow-up observations of 
these candidates are well worth exploring, especially for the single- 
transit candidates, as mass limits could be derived for most of them, 
allowing us to distinguish their real nature. This is particularly im- 
portant as the periods of these objects are all 5 days or less. If brown 
dwarfs are confirmed among the sample, then they would fall into 
the so-called "brown dwarf desert". This in turn might indicate that 
the low-mass object population in this field differs from that of the 
solar neighbourhood. 

Such a result would be interesting given the lack of transiting 
planets (and brown dwarfs) found in the old, metal-poor globular 
cluster, 47 Tuc. lBrown et alJ 1200 lh concluded that the absence of 
planets might be explained by the low metallicity, and/or crowded 
environment serving to disrupt planetary formation and evolution. 
NGC 6819 is comparatively metal-rich and provides a different en- 
vironment in which to study the importance of these factors. Our 
rough estimate suggests that we should have detected about 1 1 tran- 
siting planets in these data, if hot Jupiters are as common as they are 
in the solar neighbourhood. Of course, the transit method favours 
stars of later spectral type than the RV technique, so it is possible 
that planetary frequency decreases for later spectral type. We will 
discuss the significance of this result in detail in an forthcoming 
paper. 
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